栅格二值化后获取分布质心

# 一种是采用数据质心分布的方法：直接进行数据的二分类，这样所有的数据值都是一致性的，可以使用空间坐标系求均值的方式；
# 另外一种方式是采用加权求均值的方法，考虑到数据

解决shp文件中自交报错的问题：

library(rgeos)
sp.r1 <- rasterToPolygons(raster(xmn = -125, xmx = -70,
                                ymn = 18, ymx = 50,
                                resolution = 3))
## 下面这一部分是必须的：
data(wrld_simpl)
polygone1 <- gBuffer(sp.r1, byid=TRUE, width=0)
polygone2 <- gBuffer(wrld_simpl, byid=TRUE, width=0)
clip2 <- rgeos::gIntersection(polygone1, polygone2, byid=TRUE, drop_lower_td = TRUE)